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ABSTRACT 

We demonstrate that for several of the gravitational lens models used to describe 
galaxies, there exists a quantity we dub the magnification invariant, equaling the sum of 
the signed magnifications of the images, that is a constant when the image multiplicity 
is a maximum. This invariant is independent of most of the model parameters, and is 
independent of the source position as long as the source lies inside the caustic. It is not 
necessary to solve the lens equations to compute this invariant. For quad lenses, it is 
usually easy to assign image parities and thus one can compute the sum of the signed 
fluxes, compare with the magnification invariant of different models, and infer the model 
magnification factors of the images without fitting. We evaluate the invariant for simple 
galaxy models, apply it to known cases of quadruple lenses, and discuss implications for 
the ability of these lens models to reproduce observed image brightnesses. 

Subject headings: gravitational lensing-galaxy models-individual(CLASS 1608+656, 
H1413+117). 
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1. INTRODUCTION 

Gravitational lensing is a well-known phenomenon 
that has emerged as a useful tool in modern astro- 
physics. Two important classes of gravitational lenses 
are galaxy lenses, which can produce multiple images 
of the background source with typical image separa- 
tions of arcseconds (Keeton and Kochanek, 1996), and 
lensing by compact objects in our own and nearby 
galaxies (eg Alcock et al., 1997 and Alcock et al., 
1998). The latter class, termed "microlenses" since 
typical image separations are too small to permit indi- 
vidual image resolution, are modeled as point masses. 
Models of galaxy lenses can be much more compli- 
cated, since the mass distributions of galaxies, al- 
though not well understood, are certainly far more 
complex than those of the starlike microlenses. Since 
galaxies are thought to be dominated by dark halos 
in roughly spherical, isothermal distributions, the first 
and most important part of a galaxy lens model is an 
isothermal term. Variations on this can add elliptic- 
ity and/or shear into the mass distribution, or into 
the lens potential itself (Kochanek, 1991; Kassiola & 
Kovner, 1993; Kormann et al., 1994; Keeton et al., 
1997; Witt and Mao, 1997; Keeton and Kochanek, 
1998). The latter set of models are generally much 
easier to deal with analytically, so we will concentrate 
on those. 



Witt & Mao (1995) showed that for a binary mi- 
crolensing system, when the image multiplicity is a 
maximum (five in this case), then the sum of the 
signed magnifications of all 5 images is always unity. 
The sign of the magnification of an image is merely 
the parity of the image, +1 for ordinary images and 
-1 for inverted images, (see Schneider, Ehlers & Falco 
(SEF), 1992 for a thorough review of lens theory and 
terminology.) Writing fj,i as the signed magnification 
of the i th image, let us define / = X^Mi- This is a 
different quantity from the total magnification, which 
is given by \fn\. Witt & Mao showed that 1 = 1, 
independent of quantities such as the Einstein radii 
of the lenses or the position of the source, as long 
as the source was inside any of the caustics. Witt & 
Mao ingeniously applied resultant theory to the lens 
equations to obtain this result. More recently, Rhie 



(1997) used a coordinate transformation to derive the 
same result, 1 = 1, much more easily, and to show 
that 1 = 1 for N point masses in a plane acting as a 
gravitational lens, for all nonnegative N. Thus, for mi- 
crolensing, the quantity / is an invariant, independent 



of lens model, so we will refer to it as the "magnifi- 
cation invariant" . The invariance of the sum of the 
signed magnifications of the maximum number of im- 
ages is useful for microlensing because, even though 
individual images cannot be resolved, the total mag- 
nification is measurable, and conditions upon the sum 
of the signed magnification imply conditions upon the 
sum of the unsigned magnifications, as Witt & Mao 
describe. Galaxy lensing presents us with the reverse 
situation - although the total magnification cannot 
be determined observationally, individual images can 
be resolved. Since for many lenses, it is easy to assign 
image parities, one can compute the sum of the signed 
fluxes. Therefore, if for galaxy lensing there are con- 
straints upon the sum of the signed fluxes, as there are 
for microlensing, these constraints should be observa- 
tionally verifiable. It is therefore interesting to ask 
whether such constraints exist for models of galaxy 
lenses. In this letter, we consider several of the com- 
monly used models for galaxy lenses, and show that 
for many of them, although not all, such a property 
is true. We then apply this knowledge to several ob- 
served lens systems, and show how the magnification 
invariant can be used to determine properties of the 
lens, without fitting. 

2. METHOD 

The method we emplo y is straightforward, and is 
based upon that used by Witt fc Mao (1995) for the 
binary point lens. For the simple lens potentials we 
consider, we write the lens equations as polynomi- 
als in the variables x and fx, where x is the angu- 
lar position and \i is the signed magnification. We 
then apply elimination theory to eliminate the im- 
age positions and obtain a single equation satisfied 
by all of the magnifications. Elimination theory is a 
well developed field within algebraic geometry, and 
there are multiple techniques and algorithms known 
to eliminate variables from simultaneous polynomial 
equations. The most popular method is to compute 
the Grobner basis of the original set of polynomi- 
als. When lexicographic ordering of the variables 
is chosen, successive basis elements will have vari- 
ables eliminated (Cox et al., 1997). A more classi- 
cal technique is to compute the generalized resultant. 
Sylvester, in the 19th century, developed a formula to 
take the resultant of three polynomials in three vari- 
ables (Gelfand et al., 1994), which is the category into 
which the lens equations fall. Whatever technique is 
used, the result is a single polynomial equation in the 
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magnification, whose roots are the signed magnifica- 
tions of the images. Thus, the degree of the final 
polynomial equals the maximum number of images. 
We should remark that it is not always possible to 
perform such a procedure for all lens models, and it 
is important to verify that the degree of the final poly- 
nomial does equal the maximum number of images. 
Since the sum of the roots of a polynomial equals the 
ratio of the second coefficient to the leading coefficient 
(up to a sign) , to compute the sum of the magnifica- 
tions, we need only compute this polynomial. This 
ratio of coefficients then gives us the magnification 
invariant. It is not even necessary to solve the lens 
equations! This technique was first set out by Witt 
and Mao in their 1995 paper on binary point lenses. 

3. LENS MODELS 

We now apply this technique to several popular 
galaxy lens models. Galaxy lens models, as stated 
earlier, generally consist of variations on the singu- 
lar isothermal sphere (SIS). The lens potential for an 
ordinary SIS is iJj — br, where 6 is the Einstein ra- 
dius (which sets the physical scale of the problem) 
and r — \/ x 2 + y 2 is the distance to the center of the 
lens (SEF, 1992). Galaxies are not perfectly spheri- 
cal, however, and thus ellipticity and/or shear terms 
are added to the lens potential to account for this (see 
Blandford & Kochanek, 1987; Kochanek, 1991; Kas- 
siola & Kovner, 1993; Witt & Mao, 1997). For the 
models we consider, we let 7 measure the strength of 
the shear, and e be the ellipticity. We choose with- 
out loss of generality coordinates centered on the lens, 
and oriented along the ellipticity axis or direction of 
shear. The models are : 

1. SIS + elliptical : V = br + jbr cos(26») 



2. SIE : tp = bR = by/(l - e)x 2 + (1 + e)y 2 

3. SIS + external shear: tp = br + %r 2 cos(26>) 

4. SIE + external shear: ip = bR + \r 2 cos 2(8 - 7 ) 

Here, SIE stands for "singular isothermal ellipse" . 
This is not the isothermal elliptical mass distribu- 
tion described by jKassiola fc Kovner (1993)) and Kor- 



mann et al. (1994), but is instead an elliptical poten- 



tial, corresponding to the singular isothermal clli pti- 
cal potential (SIEP) of |Kassiola fc Kovner (1993)| . In 



model 4, the parameter 7 describes the orientation 
of the external shear relative to the ellipticity axes. 
We consider here only singular models since lensed 
systems generally consist of even numbers of images 



(Kochanek, 1991), but later we will discuss one non- 
singular model. 

As an example of the procedure, we calculate / 
for model 1, SIS+elliptical. For convenience, we use 
an orthonormal basis {e^e^} = {■§-, ^ n polar 

coordinates, with the source at (s,9 s ) and image at 
(r, 6), the time delay becomes 

t = i(r 2 + s 2 - 2rscos(0 - S )) - br - 7 6rcos(20) 

Stationarity of the time delay gives us the lens equa- 
tions, the solutions of which are the image positions: 

dr 

T~ = r ~ 
or 

ldr 
rdd 



T.f 



To 



6(1 + 7 cos(2<9)) = 



ssin(0-0 s )+2 7 6sin(20) = 



Note that the 9 equation is independent of r (Kassiola 
and Kovner, 1995). The components of the Hessian 
are 

' 1 
1 - b - + 3 7 £ cos(20) 

where we use semicolons to denote covariant deriva- 
tives. Thus, the Jacobian determinant for this lens 
(see appendix) is 

IIJII = 1 - - + 3 7 -cos(26»), 
r r 

and of course the signed magnification is /i = 1/|| J||, 
which is to be evaluated at each image position (i.e. 
at all of the solutions of the lens equations). 

We now write this equation, and the two lens equa- 
tions, as polynomial equations. Let u — e lB and 



jo. 



Then the lens equations become 



7 6u 4 



s„,3 



U + 2(6 — r)u 2 + szu + 7 6 = 0, 



jbu 4 



2z 



'-u-jb = 0, 



3-fbfiu 4 + 2({ir — r — b^L)u 2 + 3 7 6/i = 0. 

Eliminating r and u, and assuming 7 7^ 0, we obtain 
the desired fourth degree polynomial in /i, which is 
too long to print here. Again, this polynomial must 
be fourth degree since there are at most four images. 
We note that the leading coefficient, 04, and the cubic 
coefficient, 03, are related by 



a 4 = -a 3 



^^(36V(9-2z 4 + 9z 8 ) 7 ; 

z 4 



z 4 ' 

-s 6 z 4 - 76 86 4 s 2 z 4 7 4 + 40966 6 z 4 7 6 ) 



8^2 

,6,4 fiN 
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and so the sum of the roots of this polynomial equals 
— a 3 /a 4 = 1. This fourth degree polynomial has ei- 
ther 2 or 4 real roots. When the source is inside the 
caustic, all roots are real and thus I = Mi = 1- 
When the source is outside the caustic, 2 roots merge 
and become a complex conjugate pair. The real roots 
correspond to visible images, while the complex roots 
correspond to spurious solutions to the lens equations, 
which are not visible images (Petters, 1993). There- 
fore, when the number of real images is the maximum 
possible (4 in this case), 7=1. 

Applying this procedure to all of the galaxy lens 
models, we obtain Table 1. We should note that 
model 4 has the surprising property that when the 
shear is oriented along the cllipticity axis (i.e. 6* 7 = 0) 
and 7 = e, then the caustic shrinks to a single point, 
located behind the lens. Although there are only two 
images for this special case, the two magnifications 
sum to give the quantity listed in Table 1. 

4. DISCUSSION 

We have shown that for several models commonly 
used to describe galaxy lenses, the predicted image 
magnifications must obey the condition that their 
sum equal a simple constant, independent of source 
position (while the source lies within the caustic) and 
often independent of model parameters. This has 
some interesting consequences. To begin, we note 
that it is easy to compute the value of this magni- 
fication invariant for a quad lens system, since the 
parities of the images alternate as one goes around 
the ring of images. Now consider, for example, the 
object CLASS 1608+656 (Keeton & Kochanek, 1996; 
Myers et al., 1995), which is a promising candidate 
for measurement of the Hubble constant Ho. By in- 
spection, we see that images A and B should have 
positive parity and that C and D should have neg- 



Table 1: Magnification invariants for the models 
shown. Here, r = \/ x 2 + y 2 , = tan _1 (y/x), 



R = y/(l - e)x 2 + (1 + e)y 2 . 


Model 


1> 


E 4 Mi 


SIS 


br 


2 


1 SIS + elliptical 


br + "fbr cos(20) 


1 


2 SIE 


bR 


2 


3 SIS + external 


br + \r 2 cos(20) 


2 

1-7^ 
2 

1-7^ 


4 SIE + external 


bR+^r 2 cos 2(0 - 7 ) 



ativc parity (using the naming scheme used by Kee- 
ton & Kochanek). Relative to B, A has brightness 
2.06±0.06, C has 0.85±0.03, and D has 0.26±0.03. 
Thus, the magnification invariant for this system is 
1.95±0.07, multiplied by image B's unknown magni- 
fication factor. Now, Seitz & Schneider (1994), ap- 
plying the Raychaudhuri equation (see MTW, 1973 
or Wald, 1984), showed that minima of the time de- 
lay are never demagnificd, that is ^ > 1. Since 
B is a minimum, its magnification factor must be 
> 1. Therefore, we immediately see that model 1 
cannot possibly reproduce CLASS 1608+656, since 
to do so would require that image B have magnifi- 
cation 0.51+0.02. The other models considered here 
are not expressly forbidden, but are quite unlikely to 
succeed, since they would require image B to have a 
magnification factor barely above the minimum pos- 
sible value. Thus we see, without the effort of mod- 
cling, that none of the simple models considered here 
have hope of recovering the positions and brightnesses 
of object CLASS 1608+656. One can either ignore 
the fluxes when fitting, if their accuracy is suspect 
(which should not be the case for CLASS 1608+656, 
as fluxes were measured in the radio), or one must 
resort to more complicated models. The problem, as 
Kochanek (1991) has emphasized, is that the com- 
plexity of model is limited by the small number of 
observables, which is only 9 for a quad lens. Note 
that the simple models here already have 5 or 7 pa- 
rameters, so there are very few parameters that can be 
added before the problem becomes underconstrained. 
Of course, these remarks apply only for fitting on the 
image plane. When fitting on the source plane, the 
images are solutions to different lens equations, and 
thus the sum of their magnifications need not be any- 
thing in particular. Indeed, when images are near the 
critical curve, the sum of the magnifications can take 
on arbitrary values if the images do not have the same 
source, and so the magnification invariant would have 
little meaning. However, for purposes other than ver- 
ifying that an object is in fact gravitationally lensed, 
one ususally fits in the image plane. 

Now, CLASS 1608+656 is a particularly difficult 
object to model, but what about other lensed sys- 
tems? Consider the "Cloverleaf" , H1413+117 (Kee- 
ton & Kochanek, 1996; Yun et al., 1997). Using the 
I band fluxes listed by Keeton & Kochanek, which 
should be the most trustworthy, we find that the mag- 
nification invariant equals 0.07 times A's magnifica- 
tion factor. We shouldn't trust this very far however, 
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since the invariant computed from the R band fluxes 
would be negative! In any case, fitting any of the 
simple models to H1413+117 would require very large 
image magnification factors for a good fit, about 15 
for model 1 or about 30 or higher for the others. This 
is using the I band data; using the R or B data would 
give an even smaller invariant, and therefore require 
even greater image magnifications. Since the fluxes of 
the four images are all comparable, this means that 
all of the images would be required to lie very close 
to the critical curve, implying a tiny caustic with the 
source inside. Since the size of the caustic is deter- 
mined by the magnitude of the shear/ellipticity, con- 
straining the caustic to be small limits the range of 
these parameters significantly. Indeed, even though 
H1413+117 would seem to be an ideal object for fit- 
ting - one might go so far as to call it a "model" 
lens system - we were unable to fit any of the listed 
models satisfactorily to it, using the Xs error function 
described in Kochanek (1991), except for model 4. By 
satisfactory fit, we mean x 2 l^dof — 1- Model 4 was 
able to give a good fit due to the feature mentioned 
above, that when the external shear and ellipticity 
arel related in a specific way, the caustic shrinks to a 



point, even though the ellipticity may be significant. 
This feature enables model 4 to attain the large mag- 
nifications necessary, while possessing the ellipticity 
required to fit the image geometry. The other mod- 
els, which cannot be elliptical without having a large 
caustic, are thus doomed to failure. Our difficulty in 
fitting this object is not much of a revelation; previ- 
ous efforts at fitting H1413+117 using simple models 
such as the ones discussed here are unable to recover 

For examp le, 
at- 
do 



neither 


Kayser et al. (1990) 


nor 


Kochanek (1991) 


tempted to fit the brightnesses. 


Yun et al. (1997) 



fit the flux ratios, using a variant of model 4 with 
extra parameters. Their model is similar enough to 
model 4 to possess the same "shrinking caustic" fea- 
ture, and not surprisingly their best fit has a tiny 
(tangential) caustic barely enclosing the source QSO, 
with the associated large image magnifications of or- 
der 50. The prerequisite that the caustic be tiny for 
a good fit makes sense in terms of the magnification 
invariant. 

H1413+117 also illustrates another use of the mag- 
nification invariant. By merely adding up four num- 
bers, we immediately learn that the image magnifi- 
cations must be large. Yun et al. point out that this 
means that the curvature of the time delay surface 



at the ring of images is small, meaning that the rel- 
ative time delays between images must be small. We 
have seen, however, that the magnifications must be 
large for all models, and therefore the predicted rela- 
tive time delays between images will be small, inde- 
pendent of model. This is useful knowledge for Hq 
determination. 

We have so far discussed only singular models. 
Lens potentials can be rendered nonsingular by the 
introduction of a core radius. Altering the models 
discussed here to make them nonsingular destroys the 
invariance of the sum of the magnifications. However, 
the size of the core radius term is limited by the ab- 
sence of a central image from observed lens systems. 
For small enough core radius the sum of the magnifi- 
cations of the remaining four images will be close to 
the value for the singular model. Thus, our results 
should still apply for the cases of interest, where no 
central image is observed. For models based upon el- 
liptical mass distributions, one can see from numerical 
calculations that the sum of the magnifications of the 
maximum number of images is not a constant. How- 



ever, for typical parameter values, such as in Jackson 
et al. (19981 , [Myers et al. (1995)| or [Nair (1998)f we 



find that the sum of the magnifications of four images 
is between 2 and 3, so this model should be marginally 
better at fitting objects like CLASS 1608+656. Mod- 
els with multiple galaxies (e.g. Hogg & Blandford, 
1994) can have a much higher magnification invari- 
ant, however they have a maximum number of images 
not equal to four. 

The magnification invariant, we have seen, is easy 
to calculate for observed quad lenses, since one merely 
adds up four numbers. Doing so gives us useful infor- 
mation on the ability of lens models to fit the observa- 
tions, even before fitting. We have seen that for small 
invariants, such as H1413+117, large magnifications 
are necessary. On the other hand, invariants that are 
large, such as for CLASS 1608+656, can help rule out 
models a priori. It is interesting that two lenses re- 
cently discovered, B1933+503 (Sykes et al., 1997) and 
B0712+472 (Jackson et al., 1998) both have large in- 
variants and should, like CLASS 1608+656, be hard 
to fit using the simple models. We hope that mod- 
els b ased on realistic mass distri butions in galaxies, 
e.g. Keeton fc Kochanek (1998) , willl fare better at 
describing these objects. 
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APPENDIX 

In an arbitrary frame, it is not obvious how the 
Jacobian of the lens mapping is related to the Hes- 
sian of the time delay. We can easily relate these 
two by noting that in the usual Cartesian coordinates, 
J 1 ^ — r' l j, i.e. that the components of the Jacobian 
matrix equal the components of the Hessian. Note, 
however, that the Hessian in this equation has been 
contracted with the metric tensor, since the Hessian is 
a (^-tensor while the Jacobian, a linear transforma- 
tion taking 1-forms to 1-forms, must be a (J) -tensor. 
Of course, in orthonormal frames such as the one used 
above, or Cartesian coordinates, the metric tensor is 
merely the identity matrix. Knowing the relation be- 
tween the Jacobian and Hessian in Cartesian coordi- 
nates tells us their relationship in arbitrary frames. 
Let dx % be the 1-form dual to the basis vector q in 
the image plane, and similarly for ds % in the source 
plane. Then 

, ; ds 1 dsi dx k , t 
ds = tt-tt-e: — ~dx l 
dsi dx k dx l 

where unhatted components are Cartesian, and hat- 
ted components are in arbitrary frames. Substituting 
the Cartesian lensing Jacobian 



ds 3 
dx k 



dx 3 
g x m 



m 

A dx k 



into the previous equation gives us the lensing Jaco- 
bian expressed in terms of the Hessian computed in 
an arbitrary frame. 



J' 



ds 1 dx 3 .£ 



For orthonormal frames, the Jacobian matrices 
and are merely rotation matrices, with unit de- 

dx k J . ' 

terminant, so that the Jacobian determinant equals 
the determinant of the Hessian matrix. 
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